Semiclassical analysis of a complex quartic Hamiltonian 
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It is necessary to calculate the C operator for the non-Hermitian "PT-symmetric Hamiltonian 
H — \p 2 + ^)i 2 x 2 — \x 4 in order to demonstrate that H defines a consistent unitary theory of 
quantum mechanics. However, the C operator cannot be obtained by using perturbative methods. 
Including a small imaginary cubic term gives the Hamiltonian H = ±p + ^ /i x + igx — Ax , whose 
C operator can be obtained perturbatively. In the semiclassical limit all terms in the perturbation 
series can be calculated in closed form and the perturbation series can be summed exactly. The 
result is a closed-form expression for C having a nontrivial dependence on the dynamical variables 
x and p and on the parameter A. 
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In this paper we consider a quantum system described 
by the Hamiltonian 



H = 



\v 2 



12 2 

X 



igx 



-Xx 4 



(1) 



where g is real and nonzero and A > 0. Note that the 
potential is complex and that when A is positive the po- 
tential is unbounded below. This Hamiltonian is not Her- 
mitian in the conventional sense, where Hermitian con- 
jugation is defined as combined transpose and complex 
conjugate. Nevertheless, the eigenvalues E n are all real, 
positive, and discrete. This is because H possesses an 
unbroken VT symmetry 0, Q , which means that H and 
its eigenstates ip n {%) are invariant under space-time re- 
flection. Here, V denotes the spatial reflection p — > — p 
and x — > — x, and T denotes the time reversal p — ► —p, 
x — > x, and i — > —i. 

Many PT-symmetric quantum-mechanical Hamiltoni- 
ans have been studied in the recent literature 0, S S @ • 
However, the Hamiltonian is special because when 
g ^ the boundary conditions on the eigenfunctions may 
be imposed on the real- a; axis, as opposed to the inte- 
rior of a wedge in the complex- a; plane, as we will now 
show: The quantization condition satisfied by the eigen- 
functions requires that ijjn{x) must vanish exponentially 
in a pair of wedges in the complex- a; plane. These wedges 
are symmetrically placed with respect to the imaginary 
axis. The wedges have an angular opening of 60° and 
lie below the positive and negative real- a; axes with the 
upper edges of the wedges lying on the real axis. Using 
a WKB approximation we can determine the asymptotic 
behavior of the eigenfunctions, and we find that the ex- 
ponential decay of these wave functions is given by 



±V2A[ke j /3+gx^ /(4A)] 



(M-oo). (2) 



Thus, the wave functions are oscillatory on the positive 
and negative real- a; axes when g = 0. However, when g 
is nonzero the wave functions decay exponentially on the 
real axis as well as in the interiors of the wedges. Thus, 



taking g ^ allows us to treat real variable and 

to perform calculations on the real axis. 

Being able to treat x as real is crucial. The domain of 
the eigenfunctions of H = ip 2 + -|/i 2 a; 2 — Aa: 4 is the inte- 
rior of a pair of 60°-wedges in the lower-half a;-plane. 
Under space reflection x —> —x, this domain changes 
to the interior of a pair of 60°-wedges in the upper-half 
plane. Therefore, this Hamiltonian is not parity symmet- 
ric. However, when g ^ 0, the domain of the eigenfunc- 
tions of H in ||TJ includes the real- a; axis. Thus, on the 
real- a; axis, the parity operator V commutes with the —x A 
operator. This fact enables us to perform in this paper a 
perturbative calculation of C. The C operator is needed 
to formulate a consistent quantum theory described by 
the non-Hermitian Hamiltonian l[T)l. 

To make sense of H in JQ) we must identify the Hilbert 
space spanned by the eigenfunctions of H and then con- 
struct for this space an inner product that is positive defi- 
nite. As shown in Ref. 0], an inner product defined with 
respect to 'PT-conjugation leads to an indefinite met- 
ric of the type investigated earlier by Lee and Wick 0. 
However, an inner product defined with respect to CVT- 
conjugation leads to apositive definite metric, and hence 
positive probabilities 0, @ • Here C denotes a linear oper- 
ator analogous to the charge operator in particle physics. 
The operator C commutes with the Hamiltonian and its 
square is unity, so its eigenvalues are ±1. Because C com- 
mutes with H, the time evolution of the theory is uni- 
tary; that is, the norm of a vector is preserved in time. 
Given the operator C, we can construct the positive op- 
erator = CV, which can, in turn, be used to construct 
by a similarity transformation an equivalent Hamiltonian 
H = e~®/ 2 He®/ 2 . The Hamiltonian H is Hermitian in 
the conventional sense 0, but it is a nonlocal function 
of the operators x and p. 

Thus, the key step in formulating a consistent quan- 
tum theory based on the Hamiltonian JJJ is to calculate 
the operator C. When A = 0, one can use perturbation 
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theory to calculate the C operator as a series in pow- 
ers of g 9]. However, for the more interesting case of a 
negative quartic interaction (g = 0, A > 0), a perturba- 
tive calculation of C using conventional Poincare asymp- 
totics fails because to all orders in powers of A the oper- 
ator Q vanishes. Only nonperturbative techniques such 
as hyperasymptotics (asymptotics beyond all orders) [loj 
can be used to find the C operator for the Hamiltonian 
H = \p 2 + \p?x 2 - Ax 4 . 

The analysis in this paper is based on the observation 
that when g ^ 0, no matter how small, it is possible to use 
perturbative methods to calculate C. Our perturbative 
calculation is organized as follows: First, we introduce 
the small positive parameter e into the Hamiltonian i[TJl 
and consider 



H = \p 2 + \[i 2 x 2 + iegx 3 



l Xx 4 



(3) 



We seek a perturbation series in powers of e. The coeffi- 
cient of e" in this perturbation series is complicated, and 
thus our second step is to simplify the coefficient by mak- 
ing a semiclassical approximation in which we only retain 
leading order terms in Planck's constant h. The result is 
a series in powers of g, and since we may take g arbitrarily 
small, our third step is to simplify the coefficient further 
by omitting all contributions from higher powers of g. 
The resulting infinite series can then be summed exactly 
and in closed form. Once the summation is performed, 
our fourth step is to set e = 1 to obtain the semiclas- 
sical approximation to C for the Hamiltonian The 
Hamiltonian J3J was first considered by Banerjee In 
Ref. [Tlj the first seven terms in the perturbation expan- 
sion for the C operator are calculated (but not in the 
semiclassical regime). 

We begin our analysis by recalling that the C oper- 
ator can be expressed in the form C = e^V, where 
Q = Q(x,p) is a function of the quantum dynamical vari- 
ables x and p. In earlier work we showed that C can be 
determined by searching for an operator that satisfies the 
following three conditions fl2j |: 

(i) [C,VT] = 0, (ii) C 2 = 1, (iii) [H,C] = 0. (4) 

Substituting C = e^V into (i), we obtain the condition 
Q(x,p) = Q(—x,p), so Q(x,p) is an even function of 
x. Substituting C — e^V into (ii), we get Q(x,p) — 
—Q(—x,—p). Since Q(x,p) is even in x, it is odd in p. 
Finally, condition (iii) reads 



[H, e Q T] = 0. 



(5) 



Our objective is to determine the expression for the 
operator Q(x,p), when H is given by 0- Let us write 
the Hamiltonian @ in the form 



H = H + eH x + e 2 H 2 , 



(6) 



where Ho is the Harmonic oscillator Hamiltonian, Hi = 
igx 3 , and Hi — —Ax 4 . The commutation relation (J5J 



then implies 

H e Q V - e Q VH Q + e {H x e Q V - e Q VH l ) 

+ e 2 (H 2 e Q V - e Q VH 2 ) = 0. (7) 

Under parity we have 

VH V = H Q , VH 1 V = -H 1 , and VH 2 V = H 2 .{8) 

As noted above, H 2 for x real commutes with V because 
g is nonzero. 

Substituting these relations into J7J) and multiplying V 
from the right, we obtain 

e Q H ~H e Q = e (e Q H x + H ie Q ) 

-e 2 {eQH 2 - H 2 eQ) . (9) 

We then multiply by on the left and get 

H a -e- Q H e Q = e (H x + e~ Q H ie Q ) 

-e 2 (H 2 - e- Q H 2 eQ) . (10) 

In order to analyze (|10fl we make use of the Campbell- 
Baker-Hausdorff relation 

e- Q He Q = H+[H,Q] + &[[H,Q],Q] 

+^[[[H,Q],Q],Q} + --- (11) 

and the fact that Q can be expanded as a power series in 
c: 

Q = eQ 1 +e 3 Q 3 + e 5 Q 5 + -- - . (12) 

Substitution of (|12f) into the right side of (|11|) yields 

e- Q He Q = H + e[H,Q 1 } + %(?[[H,Q],Q] 
+e 3 ([H,Q 3 ] + i[[[ff,Qi],Qi],Qi]) 
+£ 4 (i[[[ ff 'Q3],Qi] + iT[[ff,Qi],Q3] 

+±[[[[H,Q 1 ],Q 1 ],Qi],Q 1 ]) 
+e 5 {j l [[[[[H,Q 1 },Q 1 },Q 1 },Q 1 ],Q 1 ] 
+j [ i[[H ) Q 1 ],Qx],Q3] + ^[[[H,Qi],Qa],Qi] 
+U[[H,Qs},Qi},Qi} + [H,Q 5 })+--- . (13) 

Inserting the expansion (|13JI into (|10[1 and equating co- 
efficients of powers of e, we obtain the following set of 
identities: 

[Q u H ]=2Hi, 

[Q 3 ,H ] = ^[[[ffo,Qi],Qi],Qi] + if[[#i,Qi],Qi] 
+[H 2 ,Qi], 

[Q 5 ,H ] = i[[[[[i? ,Qi],Qi],Qi],Qi],Qi] 
+ i[[[[H 1 ,Q 1 ] 1 Q 1 ],Q 1 ],Q 1 ] 
+^[[[Ho,Q 1 ],Q 1 ],Q 3 ] + ^[[[H ,Q 1 ],Q 3 ],Q 1 ] 
+±[[{H , Q 3 ],Qi],Qi] + &[[[H 3 , Qi],Qi],Qi] 
+U[Hi, Qs] , Qi] + U[Hi, Qi] , Q 3 ] + [H 2 , Q 3 ] ,(14) 
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and so on. These identities correspond to the coefficients 
of e, e 3 , and e 5 . The coefficients of the even powers of 
e are redundant because they can be derived from the 
coefficients of the lower odd powers of e. For example, 
the coefficient of e 2 gives [[H , Qi], Qi] = —2[Hi,Qi], 
which follows from the first relation in (|14|) . 

We now perform a semiclassical approximation in 
which we only retain terms to leading order in h. That 
is, we use the fact that each operator Qi in Ijl2(l has a 
semiclassical expansion of the form 



Q 



(-i) 



Qi 0) + HQ? 



ti z 



,(2) 



(15) 



and discard all but the leading terms Q\ ^ for i — 
1,3,5, . . .. Because we consider only the leading terms 
Ql , in what follows we omit the superscript and write 
Qi for simplicity of notation. 

We remark that in a semiclassical approximation, once 
a commutation relation is performed, we can regard x 
and p as classical variables and hence issues relating to 
operator ordering need not be considered. In this con- 
nection, the following relation applicable in semiclassical 
approximation is useful: 



[p a x\p c x d ] 



ih(bc - ad)p a+c - 1 x b+d - i . (16) 
This is a special case of the Poisson bracket relation 



\F( x ^G(x^\- l [ dFdG dF9G 
{F(x, P ), &{x,p)}-i [ 



(17) 



Using the semiclassical commutation relation (|lfi[) . we 
solve the first equation in (|14|l for Q\ and obtain 



Qi = 



hP 3 



12 2 
2^ P x 



(18) 



Substituting (|18fl into the second relation of l|14f> allows 
us to determine Q3 as 



4 2 A 3 



n 8 h 

4 2 g 3 



§p 5 + p 2 p 3 x 2 + \p 4 px i 



^p 5 + §/iV z 2 + 



(19) 



Similarly, substituting (fT%| and fH^I into (fT4"|) , we deduce 
that 



4 3 A 2 ff 



4 3 Ag 3 
4 3 5 5 



W + 2 mu 2 p 5 x 2 + 2p A p 3 x i + hf-L px 



fp 7 + 6fi 2 p 5 x 2 + f /iV 



| p 7 + ^2^2 + 1^3^ + ^6 



(20) 



Continuing in this manner, we can determine the per- 
turbative expansion of Q explicitly. Observe, however, 



that Qin+\ for each n = 0, 1, 2, ... is an odd polynomial 
of g of degree 2n + 1 . This follows from (|14fl if we notice 
that Hi and hence Q\ are proportional to g whereas H 
and H 2 are independent of g. Because we assume that 
the Hamiltonian has a weak cubic interaction, the 
value of the coupling g is small. Therefore, we may omit 
terms of order g 3 and higher from the expansion of Q. 
To first order in g the set of identities in (|14(l reduces to 
the following simpler set of relations: 



[H ,Qi] 

[H ,Qs] 
[Ho,Qr] 



-2H X 
[Qi,H 2 ] 
[03, W] 



(21) 



From these relations we deduce that to first order in g, 
Q7 is given by 



Qi 



4 4 A 3 <? 



\p 9 + 4 mu 2 p 7 x 2 + 6/i 4 p 5 x 4 



(22) 



and that to first order in g, Qq is given by 
4 5 A 4 ff 



Q 9 



v 20 h 



ifp 11 + 8/zVV + Wifp'x 



,2„9„2 



4„7 4 



14/i 6 p 5 a; 6 + 5^V^ 8 + ^ ll V 



8^,3^.8 1 1 ,,10,10 



(23) 



By repeating this procedure and determining Q 2 n+i for 
n = 0, 1, 2, . . ., we deduce, in general, that 



Q 



2n+l 



v 



~~ ^ 2fc (2»-fc + 2)! , fc , B _ afc+a 



E 

fc=0 



2 fc fc!(2n- 2fc + 3)! 



x p 



(24) 



To determine the semiclassical expression for Q, we 
must sum the product e 2n+1 Q 2 n+i i n n - For convenience 
we define the variables 



p?x 2 



8e 2 \p 2 



and (3 



(25) 



and write 



2„+i _ 46gp 3 ^ (2n-fc + 2)! fc 
e <y 2 n+i -^-(i 2^ U7KZ 07, , om q • ( 26 ) 



/j 4 /i 



fc=0 



fc!(2n- 2fc + 3)! 



By summing (|26|) in n and interchanging the orders of 
summation, we can express Q as 



4egp 3 
Q = ir-a 



\ - M) \ - 

^ (fc + 1)! ^ 

fe=0 V ' n=0 



fc 00 



(2n + fc + l)> „ 



fc=0 



(fc + !)!£' (2n + l)! 



/3 r ' 



OC 1 



4egp 3 

' a 4 h ^ 2n + 3' 

n— 



(27) 
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To determine the first sum in the right hand side of 
we use the identity (2n+fc+l)! = J °° dt t 2n+k+1 e"' . 
After performing the resulting integral, we get 



E 

fc=0 



(a/3)* 



(2n+k+l)\ 



2/3 ! 



1)! ^ (2n + l)\ 

1 n=a v ' 

' l-aft + Vft , 1 
In = — m ■ 



VP' 



I -aft- VP I -VP. 



(28) 



The summation on the left side of (|28|l converges to the 
right side provided that the inequality aft + VP < 1 is 
satisfied. More explicitly, this inequality reads 



P < 



2Xx z 



2eV2A 



(29) 



For e <C f, the summation on the left side of Ij28(l 
converges essentially in the entirety of the semiclassical 
phase space. An analogous conclusion follows in the limit 
A — > 0. For finite e and A, there is a parabolic region in 
the semiclassical phase space in which the operator Q 
converges. We believe that this region might be associ- 
ated with the region in which the corresponding classical 
trajectories are confined, although we have not studied 
this question. 

The second term on the right side of (|27(l gives 



E 



ft* 



^2n + i ftl \2 1-y/P 



(30) 



Note that the left side of (|30jl converges for ft < 1, which 
holds automatically if l|29|l is satisfied. 

Combining (|28|l and (|30|) and substituting (|25p. we fi- 
nally deduce that to leading order in g the semiclassical 
expression for the operator Q associated with the Hamil- 
tonian Q is 



Q 



Aegp 3 

\2(3 



1 -ml 



gu 2 V2 , M 2 - 4\e 2 x 2 - 2eV2Xp 
in 



aft + VP _ 1 
aft -VP ft 

gp 



16e 2 A§fi 



4Xe 2 x 2 + 2eV2Xp 2eXh' 



(31) 



We have performed the summation explicitly, so we may 
set e — 1 in i|31|) to obtain the corresponding result for the 
Hamiltonian (Q. This achieves our objective of finding 
the semiclassical approximation to the C operator for this 
Hamiltonian. Note that if we expand the right side of 
(131(1 for small A and then take the limit A — > 0, we recover 
l|T%)l. This is because, to first order in g, Qi is the only 
term that is not proportional to A. 

A complete analysis of the C operator for a — A.t 4 
quantum-mechanical theory would be of immense im- 
portance because it could lead to an understanding of 
its — A0 4 field-theoretic counterpart. This field theory is 
asymptotically free 0, Q, ^| and might well describe 
the Higgs sector in the standard model. Of course, the 



perturbative method used here does not apply directly to 
a pure quartic — Aa: 4 theory, which is inherently nonper- 
turbative; that is, we cannot set g = to obtain the semi- 
classical expression for Q in the — Aa; 4 theory. However, 
the work we have presented here is a first step towards 
our goal of obtaining a complete semiclassical and non- 
perturbative treatment of the — Aa; 4 theory in quantum 
mechanics. 
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